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Abstract 

In relation with experiments on correlated kaons at ^-factories, it is shown 
that the request of complete positivity is necessary in any physically con- 
sistent description of neutral kaons as open quantum systems. 



Non-standard phenomenological models incorporating loss of quantum coherence and 
entropy increase can be studied in the broad framework of open quantum systems. M-[H 
These systems can be modelled as being small subsystems in "weak" interaction with 
large environments. The reduced dynamics for the subsystem is obtained by eliminating, 
i.e. tracing over, the environment degrees of freedom; in the markovian approxima- 
tion, justified in most physical contexts, it consists of a set of completely positive linear 
transformations, with forward in time composition law (semigroup property) and entropy 
increase. These maps form a so-called quantum dynamical semigroup. 

This description is rather general and has been used to model a large variety of phys- 
ical situations, ranging from the study of quantum statistical systems, ||-[|| the analysis 
of dissipative effects in quantum optics, 0-0 the treatment of the interaction of a mi- 
crosystem with a macrosystem.0-0 It has also been recently applied to the analysis 
of the evolution and decay of the system of neutral kaons. fl0ll-[|l4J1 (For motivations and 
earlier results, see 01-[E2|). 



According to standard phenomenology, [23] kaon states can be effectively described by 



2x2 density matrices p acting on a two dimensional Hilbert space. The time evolution 
u t will transform the initial state p into the state pit) = ui t [p], at time t. A consistent 
statistical description of the initial density matrix p as a state is assured by the positivity 
of its eigenvalues that are interpreted as probabilities. Clearly, for this description to hold 
for all times, the evolution map u> t must be positive, i.e. it must preserve the positivity 
of the eigenvalues of p(t), for any t. 

Complete posit ivity[U]-|| is a more stringent condition; it guarantees the positivity 
of the eigenvalues of density matrices describing states of correlated kaons, e.g. those 
produced in 0-meson decays. States of entangled, but not dynamically interacting kaons, 



evolve according to the factorized product uj t ®^t of the single- kaon dynamical maps. [|T2f If 
oj t is not completely positive, there are instances of correlated states that develop negative 
eigenvalues; in such cases, their statistical and physical interpretation is lost. []T3J| The aim 
of this note is to show that the issue of complete positivity is not only theoretical, but can 
be given experimental relevance. Indeed, we shall explicitly find that some experimentally 
accessible kaon observables, defined to be positive, would return, in absence of complete 
positivity, negative mean values. 

The evolution in time of any kaon state p can be described in general by an equation 
of the following form: 



at 



-iH cS pit) + ip{t) Hl s + LdW)} . (1) 



The first two pieces in the r.h.s. give the standard hamiltonian contribution, while Lp is 
a linear map that encodes possible dissipative, non-standard effects. 
The effective hamiltonian, 

H cS = M- iFjl , (2) 

includes a non-hermitian part describing the kaon decay; its generic form further encodes 
possible indirect CP and CPT violations. It is convenient to introduce the even and odd 



CP eigenstates \Ki), \K 2 ) and the eigenstates of H eS 
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such that ilfefflifs) = A5I-K5) and if e fj|i^i) = Ai,|.Kl). The entries of the matrices M and 
J 1 can be expressed in terms of the two complex parameters €3, e^ and the K s , K L masses, 
m s , m L , and decay widths 75, jl, so that A5 = 7715—275/2, A^ = m L —i^ L /2. For later use, 
we introduce the following positive combinations: Ar = 75 — 7l, Am = m^ — 7715, as well 
as the complex quantities r± = V ± iAm and Ar± = Ar ± 2iAm, with T = (75 + 7l)/2. 
The form of the additional piece Lp> [p] is uniquely fixed by the physical requirements 
that the complete time-evolution u t needs to satisfy; as already mentioned in the intro- 
ductory remarks, these are semigroup property, entropy increase and complete positivity. 
Expressing the kaon state p as the 4- vector (p , Pi, P2, P3) along the basis consisting of the 
identity matrix o"o and the Pauli matrices <Tj, i — 1, 2, 3, the additional linear piece Lp 
acts as the 4x4 matrix: fO 
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where a, b, c, a, /3 and 7 are six real parameters, with a, a and 7 non negative. Further, 
these parameters must satisfy the following inequalities 



2R = a + -f-a>0, RS > b 2 
2S = a + 7-ct>0, RT>c 2 
2T = a + a-7>0, ST > (3 2 

RST >2bc(3 + R(3 2 + Sc 2 + Tb 2 , 



(5) 
(6) 
(7) 
(8) 



which are necessary and sufficient conditions for the complete positivity of the time- 
evolution map cu t generated by ([[]). 

It should be stressed that in absence of the piece Ld[p], pure states (i.e. states of 
the form \tfj)(ip\) are transformed by the evolution equation (|1|) back into pure states, 
even though probability is not conserved, a direct consequence of the presence of a non- 
hermitian part in the effective hamiltonian H e g. Only when the extra piece Lj^[p] is 
also present, p(t) becomes less ordered in time due to a mixing-enhancing mechanism: it 
produces dissipation and possible loss of quantum coherence.^ 

A rough dimensional estimate gives the magnitude of the non-standard parameters in 
(fD; they are at most of the order m 2 K /Mp ~ 10~ 19 GeV, with mx the kaon mass and Mp 
the Planck mass;[|16j, pl| , [ID , [|T4|j the magnitude of this estimate is not far from the 



Although in principle these phenomena could also contribute to the decay process of the X-mesons, 
their effects can be estimated to be negligible for any practical considerations; H the description of the 
if-meson "self-dynamics" by means of the effective hamiltonian (0) is therefore appropriate. 



value of Ar e$ and Ar e^. This allows a perturbative treatment of the evolution equation 
in (|T]). It is thus possible to compare the predictions deduced from the solutions p(t) 
of (|I|) with the available data of various observables accessible in kaon experiments. | jTT| 
Unfortunately, the accuracy of the present data, essentially from fixed target experiments, 
are not sufficient to draw a definite conclusion whether the constraints ©-(§) are obeyed 
or not. More precise results should come from the experiments on correlated kaons in the 



so-called ^-factories. 24 



Since the 0-meson has spin 1, the two neutral spinless kaons produced in a 0-decay 
fly apart with opposite momenta in the meson rest-frame; their corresponding density 
matrix pa is antisymmetric in the spatial labels. By means of the projectors onto the CP 
eigenstates, 

P 1 = \K 1 )(K 1 \, P 2 = \K 2 )(K 2 \, (9) 

and of the off-diagonal operators, 

P 3 = \K 1 )(K 2 \, P 4 = \K 2 )(K 3 \, (10) 

we can write 

PA = \{Pl ®P2 + P2®Pl-Pz®Pi-P±® P 3 ) ■ (11) 

The time evolution of a system of two correlated neutral if-mesons, initially described 
by Pa, can be analyzed using the sigle if-meson dynamics so-far discussed. Indeed, as in 
ordinary quantum mechanics, it is natural to assume that, once produced in a decay, 
the kaons evolve in time each according to the completely positive map u t generated 
by (p. As already remarked, this assures that the resulting total evolution map is still 
positive for any time and of semigroup type. We stress that this choice is the only natural 
possibility if one requires that after tracing over the degrees of freedom of one particle, 
the resulting dynamics for the remaining one be completely positive, of semigroup type 
and independent from the initial state of the other particle. 

Within this framework, the density matrix that describes a situation in which the first 
if -meson has evolved up to proper time t\ and the second up to proper time t 2 is given 
by: 

PA(tuh) = (wt! ®uj t2 )[p A ] 

= ^[P 1 (t 1 )®P 2 (t 2 ) + P 2 {t 1 )®P 1 {t 2 ) - Psih) ® Pt(t 2 ) - Ptih) ® P 3 (t 2 j\ (12) 

where Pi{t\) and Pi(t 2 ), i — 1,2, 3, 4, represent the evolution according to (|1|) of the initial 
operators @, (0), up to the time t\ and t 2 , respectively. In the following, we shall set 
ti—t 2 — t, and simply call pAif) = PA(t,t) the evolution of flTT| ) U P to proper time t. 

The statistical description of Pa(£) allows us to give a meaningful interpretation of the 
quantities 

V l] {t) = r Ti[p A {t)P®P 3 ] , *,J = 1,2, (13) 

as the probabilities to have one kaon in the state \Ki) at proper time t, while the other 
is in the state \Kj) at the same proper time. When i,j = 3,4, the quantities Vijit) are 



complex and do not represent directly joint probabilities. However, as we shall see, they 
can still be obtained from data accessible to experiments. 



Following the discussion in |L2| , up to first order in all small parameters, one finds: 
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(14) 
(15) 
(16) 
(17) 
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(19) 
(20) 



The remaining quantities Vij(t) can be derived from the previous expressions by using 
the following properties: 



V l3 (t) = P* 4 (t) . 



i,j = 1,2,3,4 
2 = 1,2, 



(21) 
(22) 
(23) 



Putting a = 6 = c = a = /3 = 7 = 0, one obtains the standard quantum mechanical 
effective description that evidentiates the singlet-like anti-correlation in pA(t): Vu(t) = 0. 

We emphasize that none of the above expressions contain the standard CP, CPT- 
violating parameters eg, e^. This fact makes possible, at least in principle, a direct deter- 
mination of the non-standard parameters irrespectively of the values of eg, e^; one needs 
to fit the previous expressions of the quantities V%j (t) with actual data from experiments 
at 0-factories. 

To be more specific, we shall now explicitly show how the quantities V%j can be directly 
related to frequency countings of decay events. First, notice that, given any single-kaon 
time-evolution p \— > p(t), the matrix elements of the state p(t) at time t can be measured 
by identifying appropriate orthogonal bases in the two-dimensional single kaon Hilbert 
space. The choice of the CP-eigenstates \K\) , \K 2 ) is rather suited to experimental 
tests. Indeed, since a two-pion state has the same CP eigenvalue of \Ki), the probability 
Pt{K\) = (Ki\p(i)\Ki) of having a kaon state K\ at time t is directly related to the 
frequency of two-pion decays at time t. Possible direct CP violating effects, the only ones 
allowing K2 — > 2ir, can be safely neglected; they are proportional to the phenomenological 
parameter e', that is expected to be very small, \e'\ < 1O~ 6 .0 

2 The recent measures of TZe(e' /e), pqj,|26| seem to confirm these theoretical expectations. [B7J 



On the other hand, while the decay state 7r 7r°7r° has CP = — 1, the state 7r + 7r~7r° 
may have CP = ±1. Thus, the probability P t (K 2 ) = (K 2 \p(t)\K 2 ) to have a kaon state 
K 2 at time t is not as conveniently measured by counting the frequency of the three-pion 
decays. To avoid the difficulty, the following strangeness eigenstates can be used: 

= uf.) + \ Ka ) m = \ Kl} -m (24) 

Then, the probabilities P t (K°) = (K°\p(t)\K^)_ and P t (K°) = (K°\p(t)\K°), that the 
kaon state at time t be a K°, respectively a K°, may be experimentally determined by 
counting the semileptonic decays K \— > ti~C + u, respectively K° \— > ti + C~V, the exchanged 
decays being forbidden by the AS" = AQ rule. (In the Standard Model, this selection 
rule is expected to be valid up to order 10~ 14 .p8|) Further, the probability P t (K 2 ) = 
(K 2 \p(t)\K 2 ) of having a kaon state K 2 at proper time t can be expressed as P t (K 2 ) = 
P t (K°) + PtiK 5 ) - P t (I<i), by writing 

\K 2 )(K 2 \ = \K°)(K°\ + \W)(W\ - |#i>(#i| . (25) 

Hence, Pt(K 2 ) can be measured by counting the frequencies of semileptonic decays and 
of decays into two pions. 

In order to measure the off-diagonal elements {Ki\p\K 2 ), (K 2 \p\Ki), we use the identity 

\K°}(K°\ - \K^)(K^\ = \Ki)(K 2 \ + \K 2 ){K X \ , (26) 

and extract |i^i)(i^2| from it. To do this, we need a third orthonormal basis of vectors 
whose projectors are measurable observables in actual experiments. An interesting possi- 
bility is based on the phenomenon of kaon- regeneration (see p9f and references therein). 



The idea is to insert a slab of material across the neutral kaons path; the interactions of 

the K , K° mesons with the nuclei of the material "rotate" in a known way the initial kaon 

states entering the regenerator into new ones. As initial states, consider the orthogonal 

vectors 

~ \K x )-rf\K*) ~ y_m + \K 2 ) 
\ Ks > = / , = ' \ K l) = / , = . I 27 ) 

Vi + N 2 vi + l^/l 2 

where rj is a complex parameter which depends on the regenerating material. By carefully 
choosing the material and the tickness of the slab, one can tune the modulus and phase of 
i] in such a way to completely suppress the Kl component and to regenerate the Ks state 
into a Ki, just outside the material slab. Thus, the probability P t (Ks) = (Ks\p(t)\Ks) 
that a kaon, impinging on a slab of regenerating material in a state p(t) at time t, be 
a Ks, can be measured by^ counting the decays into 2tc just beyond the slab. Now, the 
projector onto the state \K$) reads 



1 ir^wr^l , \V 
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11 ^Ki){K 2 \--\-\K 2 ){Ki\. (28) 



Then, from (H)-© and (11) it follows that 



\K X ){K 2 \ = f 1 \K 1 )(K 1 \ + f 2 \Ks)(Ks\ 

+ h\K°)(K°\ +ft\K°)(K°\, (29) 



where 



1 — \r]\ 2 1 + \rj\ 2 

2iTm{r\) ' 2ilm(r]) 

! I 2, ^ I I z i =fc 

M -ry = |g| + ?? 

liXvnij]) ' 2iXm{rj) 

In this way, the determination of the off-diagonal elements of p(t) amounts to counting 
the frequencies of decays into two pions with or without regeneration and the frequencies 
of semileptonic decays: 

(KMt)\K 2 ) = fiPt(K 1 ) + f 2 Pt(Ks) 

+ f 3 P t (K°) + UP t (K°). (32) 

The application of these results to the case of correlated kaons is now straightforward. 
For sake of compactness, we identify the various kaon states with the projections Q^, 
ji = 1, 2, 3,4, where: 

Q 1 = \K } ){Ki\, Q 3 = \K°}(K°\, (33) 

Q 2 = \K S )(K 3 \, Q,= \K^){W\. (34) 

As discussed, these operators can be measured by identifying 2n final states, in absence 
and presence of a regenerator (Q% and Q 2 ), and semileptonic decays (Qz and Q 4 ); the 
same holds for the projectors in (^) and (|TUD, since: 

(35) 
(36) 

(37) 

Further, we denote by P t (Q^,Q u ) the probability that, at proper time t after a 0-decay, 
the two kaons be in the states identified by Q^ and Q u , respectively. Then, the determi- 
nation of the quantities Vij(t) reduces to measuring joint probabilities, i.e. to counting 
frequencies of events of certain specificied types. Indeed, one explicitly finds: 

V u (t) = Pt(Qi,Qi) , (38) 

V 12 (t) = P t (Q u Q 3 ) + P t (Q u Q 4 ) - P t (Q u QO , (39) 

Vi3(t) = J2U p t(Qi,Q»), (40) 

6 
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4 
= Pa = Z^ Jv ^»m • 



v 22 {t) = p t (Q u go + p t (Q 3 , g 3 ) + p t {Q 4 , g 4 ) 



+2 



Pt(Qs, Qa) - P t (Qi, Q 4 ) - Pt(Qi, Qa) 



p 23 (t) = £ /„ [Pt(Q 3 , Qm) + ^t(g 4 , g M ) - p*(Qi, g, 

/U=l i/ = l 



(41) 
(42) 

(43) 



As a result of the previous analysis, the inconsistencies of models without complete 
positivity, besides being theoretically unsustainable, turn out to be experimentally expos- 
able. Let P u and P v project onto the correlated states 

\u) = -±=(\K 1 )®\K 1 ) + \K 2 )®\K 2 )), (44) 

\V) = -^(l^!)®!^) + 1^)®!^)) . (45) 

The averages of these two positive observables with respect to the state pA(t) read 



U{t) = Tr[p A (t) P u ] = (u\p A (t)\u) 

= ^(Pn(t) + V 22 {t) + V 33 (t) + Vu{t)) 

V(t) = Tr[p A (t)P v ] = (v\p A (t)\v) 

= - (V 12 (t) + V 21 (t) + P 3 4 (t) + P43 (t) ) 



(46) 



(47) 



and, as explained before, can be directly obtained by measuring joint probabilities in 
experiments at ^-factories. On the other hand, (|I"4T)-(|2"0"1) give, up to first order in the 
small parameters, 



U(t) = e 



-2Ft 



-^ sinh(tAr) + —- (l - cos(2tAm) 
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V(t) = e- 2Tt (a + a--f)t . 
Thus, W(0) = V(0) = 0, whereas 



dU(0) 
dt 



a + 7 — a 



dV(0) 
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a + a — 7 



(48) 
(49) 

(50) 



are both positive because of conditions (^) and (|7|). More in general, the mean values 
( f46|) and ([|7]) are surely positive, for the complete positivity of the single-kaon dynamical 
maps uj t implies p(t) = J2j Vj(t) P^j (^);[Q] _ where the Vj(t) are 2x2 matrices such that 



J2j Vj(t)Vj(t) is a bounded 2x2 matrix.^ Then, the complete evolution p A —> Pa(<0 = 
J2i,jWi(t) ® Vj(t)] Pa Wi(t) ® Vj r(t)] will never develop negative eigenvalues. 

On the other hand, if the single-kaon dynamical map u t is not completely positive, 
inconsistencies may emerge. As an example, take the phenomenological models studied in 
[ PHI , |22| , where the non-standard parameters a, b, c are set to zero and a ^ 7, 07 > (3 2 . 
The corresponding dynamics is not completely positive: the inequalities @-@ are in 
fact violated. In this case, one still has U(0) = V(0) = 0, but 

Therefore, one of the mean values fl46|), fl4~7|) starts assuming negative values as soon as 
t > 0. The inconsistence is avoided only if a = 7, which is a necessary condition for 
getting back the property of complete positivity. We stress that planned set-ups at 0- 
factories can measure the two mean values in (|46D and (f47D and therefore clarify also 
from the experimental point of view the request of complete positivity for the dynamics 
of neutral kaons. 
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